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The isoscalar giant monopole resonance (ISGMR) in Cd, Sn and Pb isotopes has been studied 
within the self-consistent Skyrme Hartree-Fock+BCS and quasi-particle random phase approxima- 
tion (QRPA). Three Skyrme parameter sets are used in the calculations, i.e., SLy5, SkM* and SkP, 
since they are characterized by different values of the compression modulus in symmetric nuclear 
matter, namely -R'oo=230, 217, and 202 MeV, respectively. We also investigate the effect of different 
types of pairing forces on the ISGMR in Cd, Sn and Pb isotopes. The calculated peak energies and 
the strength distributions of ISGMR are compared with available experimental data. We find that 
SkP fails completely to describe the ISGMR strength distribution for all isotopes due to its low 
value of the nuclear matter incompressibility, namely K 00=202 MeV. On the other hand, the SLy5 
parameter set, supplemented by an appropriate pairing interaction, gives a reasonable description 
of the ISGMR in Cd and Pb isotopes. A better description of ISGMR in Sn isotopes is achieved by 
the SkM* interaction, that has a somewhat softer value of the nuclear incompressibility. 
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I. INTRODUCTION 

The compression modulus associated with the nuclear 
many-body systems plays an important role in the de- 
scription of the structure of finite nuclei, in the dynam- 
ics of heavy-ion collisions, and in the physics of neu- 
tron stars and core-collapse supernovae Q . For more 
than thirty years, much effort has been done to deduce 
the value of the nuclear matter incompressibility both 
theoretically and experimentally. The measurements 
of the compression modes, such as the isoscalar giant 
monopole resonance (ISGMR) in finite heavy nuclei, 
have been the best tool so far to determine the value of 
the nuclear matter incompressibility Kqo- The analysis 
has mostly, but not only, been based on the distribu- 
tion of the ISGMR strength in ^°^Ph. The main results 
have been reviewed in Ref . Q . The extracted values of 
Koo are somewhat model-dependent but the accepted 
value from ^"Spb is 240±20 MeV. In general, mean field 
models, either non-relativistic or relativistic, have been 
used to extract the value of the nuclear matter incom- 
pressibility. Widely used non-relativistic Skyrme mod- 
els give a value for the nuclear matter incompressibility 
around 230 MeV. It had been claimed that relativistic 
mean field models are characterized by larger values of 
the nuclear matter incompressibility (around 250 MeV). 
However, using new fitted Skyrme forces with different 
density dependence, the authors of Ref. [1] have pointed 
out that forces with Koo = 250 MeV can also reproduce 
the ISGMR experimental data of ^osp]-, ygj.y ^^^n g^,^ 
the residual model dependence in the extracted value 
of Koo is attributed to the fact that the distribution of 
ISGMR in ^"*Pb is also sensitive to the density depen- 
dence of the symmetry energy 4- 6] . 

Recently, the distribution of ISGMR strength in Cd, 



Sn and Pb isotopes has been measured at the Research 
Center for Nuclear Physics (RCNP) at Osaka Univer- 
sity urn. These data raise a further question on the 
nuclear matter incompressibility. Namely, the Skyrme 
effective interactions with Koo ^ 230 MeV and the rel- 
ativistic mean field (RMF) models having Koo ^ 270 
MeV can reproduce the experimental ISGMR distribu- 
tion in ^"^Zr and ^osp^ very well. However, the same 
models overestimate the centroid energies of the ISGMR 
in Sn isotopes [1, [TTI - [l3t . This discrepancy shows the 
fact that the observed ISGMR in Sn isotopes is softer 
than those in ®°Zr and ^°^Pb, and this might be related 
to our incomplete understanding of surface, asymmetry 
and pairing contributions to the incompressibility of fi- 
nite nuclei. 

Skyrme models have been used in Ref. [U to in- 
vestigate the correlation between the asymmetry con- 
tribution to the incompressibility Kr and Kac, and to 
validate the extraction of Kr from the Sn data. Models 
having K^- and Koo compatible with the Sn data, and 
built within the RMF framework, arc shown to signif- 
icantly underestimate the the distribution of strength 
in 208pb By calculating the ISGMR in nuclei 

with large neutron excess, the authors of Ref. [lB| con- 
cluded that the incompressibility of neutron-rich mat- 
ter is still an important open problem. It has also been 
pointed out that superfluidity may have a sizable ef- 
fect on the incompressibility of nuclear matter and fi- 
nite nuclei [13]. This conclusion has been drawn 
by exploiting constrained Hartree-Fock or Hartree-Fock 
Bogoliubov calculations [1, [l^, [2^ in order to determine 
the inverse energy- weighted sum rule m_i, and by using 
this together with the energy-weighted sum rule toi to 
define the ISGMR centroid as i?isGMR = v^mi/m_i. It 
should be stressed that the effect of pairing on the IS- 
GMR, within the self-consistent quasi-particle Random 
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Phase Approximation (QRPA) on top of HFB, has been 
highhghted in Ref. ill|J : in this work it has been shown 
that the inclusion of pairing reduces, for some extent, 
the discrepancy between the values of K^o extracted 
from ^'^^Pb and Sn isotopes data, respectively. 

A different type of calculation was performed in Ref. 
[21I] to describe the ISGMR strength distribution in Sn 
isotopes and ^°*Pb. The theoretical models used are 
the QRPA and the quasi-particle time blocking approx- 
imation (QTBA), that includes quasi-particle-phonon 
coupling. Also in this satisfactory description of 

^•^^Pb and Sn isotopes at the same time has not been 
achieved. This calculation is not a fully self-consistent 
one. 

In this work, we employ the Skyrme QRPA approach 
on top of HF-Bardeen-Cooper-Schrieffer (HF-BCS), to 
study the ISGMR strength distribution in Cd, Sn and 
Pb isotopes. This method allows systematic and fully 
self-consistent calculations for ISGMR. All terms of the 
interaction, including the one-body and two-body spin- 
orbit and Coulomb parts, are included when the ground- 
state mean-field and the residual interaction are evalu- 
ated. Three Skyrme parameter sets, SLy5 '22], SkM* 
[23| . and SkP [2J] are used in the calculations. These 
Skyrme interactions display different values for the nu- 
clear matter incompressibility. We also compare the ef- 
fect of volume, surface, and mixed pairing interactions 
on the ISGMR properties in Cd, Sn and Pb isotopes. 
This paper is organized as follows. We will briefly report 
the main features of our Skyrme HF-BCS plus QRPA 
model in Sec. II. The results for ISGMR in Cd, Sn and 
Pb isotopes are discussed and compared with available 
experimental data in Sec. III. Section IV is devoted to 
the summary and discusses the perspectives for future 
work. 



II. THEORETICAL MODELS 

The Skyrme interaction is quite successful in the de- 
scription of nuclear properties both of ground states and 
excited states. As mentioned in the previous Section, we 
will use the QRPA approach on top of HF-BCS for our 
theoretical investigation. 

Firstly we solve the Skyrme HF-BCS equation for the 
ground state in coordinate space. The radial mesh on 
which the equations are solved extends up to 18 fm, 
and the mesh size is 0.1 fm. For all the nuclei under 
study, this radial mesh is large enough so that the results 
are stable. The pairing correlations are generated by a 
density-dependent zero-range force. 



Vpair{ri,V2) = Vb 
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where p(r) is the particle density in coordinate space 
and pq—Q.IQ ivoT^ is the nuclear saturation density. The 
value of rj is taken as 0, 0.5, or 1 for the volume, mixed, 
or surface pairing interactions, respectively. The pairing 
window, that is, the states that are taken into account 
for the solution of the BCS equations, includes five un- 
occupied orbitals above the last occupied level in the 
HF approximation. In this space, the pairing strength 



Vq is fixed by fitting the pairing gap extracted from ex- 
perimental data of odd-even mass difference by using 
the five point formula. The values of Vq that we have 
obtained are displayed in Table I. They reproduce well 
the gap associated with one typical nucleus for each iso- 
tope chain; in particular, we have chosen to reproduce 
the empirical pairing gaps of ^^^Cd (A„=1.334 MeV), 
i20Sn (A„ = 1.321 MeV), and 204pb (A„=0.841 MeV). 
For Sn and Pb isotopes, there is only neutron pairing 
because these nuclei have closed proton shells associ- 
ated with Z = 50 and 82, respectively. For Cd isotopes, 
the proton pairing also exists in principle because the 
proton number 48 is not a magic number. However, we 
have found numerically that the effect of proton pairing 
in Cd isotopes on the ISGMR strength distribution is 
very small (of the order of tenths of keV). So in our cal- 
culation we use the filling approximation for the proton 
1.99/2 state in Cd isotopes. The pairing strengths Vq 
obtained by fitting the empirical gaps are different for 
each Skyrme parameter set and for each pairing model. 

We now provide few details about the QRPA calcu- 
lations. The single-particle continuum is discretized by 
setting the nuclei in a spherical box of radius equal to 
18 fm. For every value of the quantum numbers (/, j) 
associated with the single-particle states, we include in 
the QRPA model space the unoccupied states up to the 
maximum number of nodes given by Umax = '^iast+12, 
where niast is the number of nodes of the last occupied 
state with a given {I, j). The convergence of the cal- 
culated results is checked by looking at the results for 
the energy and the strength of the ISGMR. The QRPA 
matrix equation having good angular momentum and 
parity J'^ is given by 





(2) 



where tuun is the energy of the n-th QRPA state and 
X", Y" are the corresponding forward and backward 
amplitudes, respectively. The explicit forms of the ma- 
trices A and B are given elsewhere _25ri28i] . The p- 
h matrix elements are derived from the Skyrme en- 
ergy density functional including all the terms such as 
the two-body spin-orbit and two-body Coulomb inter- 
actions. We should mention that several previous works 
devoted to the study of the ISGMR strength in Sn iso- 
topes are not fully self-consistent [III IT^. l2l| since the 
spin-orbit interaction is not taken into account in the 
residual interaction. In Ref. [2^, the authors have dis- 
cussed the self-consistency violation in HF-RPA calcula- 
tions for nuclear giant resonances: they have shown, for 
example, that the two-body spin-orbit interaction gives 
a slight repulsive contribution to the ISGMR strength 
in light nuclei, whereas, for medium and heavy nuclei, 
it produces an attractive effect on the ISGMR strength, 
so that the centroid energies are pushed downward by 
about 0.6 MeV. 

After solving the QRPA equations, various moments 
of the strength distributions can be obtained by means 
of the equation 



J E^S{E)dE, 
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lAoL/rj i: ine pairing 
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various types of pairing interactions defined 


m tq. (yj), 


in units of MeV fm'^ . For 


details see the text. 
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where S{E) = J2n \{^\F\n)\^S{E - E,,) is the strength 
function associated with the monopole operator 

F-T.'f, (4) 

i 

The constrained energy Econ, the centroid energy Ecen, 
and the scahng energy Es of the resonance are then 
defined as 




respectively. 

III. ISGMR RESULTS AND DISCUSSION 

In order to check how much the pairing correlations 
affect the ISGMR strength distributions, we show in 
Fig. 1 the ISGMR strength distributions in "°Cd and 
^^°Sn calculated by using different pairing models. The 
effective force SkM* is adopted in the particle-hole chan- 
nel. For the pairing channel, we take the surface, mixed 
and volume pairing interactions and compare the re- 
sults with those obtained within RPA (no pairing) and 
the filling approximation. We find that the surface pair- 
ing and the filling approximation give almost identical 
results. On the other hand, the volume and the mixed 
pairing predict slightly lower peak energies although the 
difference is rather small. The ISGMR centroid energy 
in "OCd (i20Sn) are 16.38 (15.80), 16.33 (15.86), 16.10 
(15.69), 16.18 (15.76) MeV, when calculated by employ- 
ing either the filling approximation, or the surface, or 
the volume, or the mixed pairing interactions, respec- 
tively. The maximum difference due to different pairing 
models is 280 (110) keV in ""Cd (^^OSn). 

A. Cd Isotopes 

Very recently, measurements of the ISGMR strength 
distributions in Cd isotopes were performed at RCNP, 
Osaka University 0. In keeping with the fact that it 
is difficult to reproduce equally well Pb and Sn iso- 
topes with a unique Skyrme force, we would like to see 
what are the results that these forces provide for the 



strength distributions in Cd isotopes. To address this 
question, we have performed the calculations for the IS- 
GMR strength in the Cd isotopes with the three differ- 
ent aforementioned Skyrme interactions together with 
various pairing models. Figure 2 displays the QRPA re- 
sults for the ISGMR strength distribution in ^^^'^^^Cd, 
calculated by using the SkP (dotted line), SkM* (dashed 
line), and SLy5 (solid line) interactions, respectively. 
The pairing force adopted in Fig. 2 is the mixed pair- 
ing interaction. For all nuclei, the strength distributions 
are concentrated in a single peak around 16 MeV, which 
exhausts almost all the energy weighted sum rule. How- 
ever, the location of the peak found with each Skyrme 
interaction is slightly different. The SkP interaction 
predicts lowest peaks while the SLy5 interaction gives 
peaks at the highest energies. The peaks obtained by 
using the SkM* interaction stays in the middle between 
the other two cases. As it is known from previous stud- 
ies, the relative position of the peaks is governed by the 
nuclear matter incompressibility associated with each 
effective interaction. 

In Fig. 3 we show the theoretical ISGMR centroid en- 
ergies in obtained by the QRPA calculations. 
The results in Figs. 3(a), 3(b) and 3(c) are obtained 
by using the SLy5, SkM*, and SkP Skyrme interac- 
tions together with various pairing interactions. Figure 
3(a) shows the ISGMR centroid energies obtained by the 
SLy5 interaction. With either the filling approximation 
or the surface pairing interaction, the calculated results 
are about 300 keV systematically larger than the exper- 
imental data. On the other hand, the calculations with 
volume and mixed pairing interaction reproduce the ex- 
perimental centroid energies within 150 keV. For the 
case of the SkM* interaction whose results are shown in 
Fig. 3(b), the conclusion is different from the one ob- 
tained by the SLy5 interaction, due to the lower value 
of the nuclear matter incompressibility (we remind that 
Koo = 217 MeV for SKM*, and that K^o = 230 MeV for 
SLy5). It seems that the filling approximation and the 
surface pairing interaction provide a good reproduction 
of the experimental centroid energies, whereas the vol- 
ume and mixed pairing interactions somewhat underes- 
timate the experimental data. Finally, it should be ex- 
pected that the results obtained with the SkP parame- 
ter set underestimate the experimental ISGMR centroid 
energies in ios-hoq^j Jj,^ pjg_ 3(c)]^ due to the lower 
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FIG. 1: The ISGMR strength distribution in (a) "°Cd and (b) ^^"Sn, calculated by using either the filling approximation, 
or the volume, or the surface, or the mixed pairing forces. The SkM* force is adopted in the p-h channel. 



TABLE II: The calculated ISGMR constrained energies (Econ), centroid energies (i?cen), and scaling energies (Es) in even- 
even ios-hgq^ isotopes are compared with the experimental data. The theoretical results are obtained in the interval 
between 10.5 and 20.5 MeV, by using the SKP, SKM*, and SLy5 parameter sets together with the mixed pairing interaction. 
The experimental data are from Ref. 0]. The values in parenthesis are the difi'erences between the theoretical values and 
the experimental data. Units are MeV. 
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SLy5 






loeCd 


16.5±0.2 


15.81(-0.69) 


16.40(-0.1) 


16.75(0.25) 






"°Cd 


16.1±0.2 


15.50(-0.60) 


16.11(0.01) 


16.46(0.36) 




y/nii/m-i 


"^Cd 


16.0±0.2 


15.35(-0.65) 


15.97(-0.03) 


16.29(0.29) 






"*Cd 


15.8±0.2 


15.21(-0.59) 


15.84(0.04) 


16.16(0.36) 
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16.3±0.1 


15.56(-0.74) 


16.18(-0.12) 


16.54(0.24) 




= mi/mo 
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15.42(-0.78) 


16.05(-0.15) 


16.39(0.19) 






"^Cd 


16.1±0.1 


15.28(-0.82) 


15.91(-0.19) 


16.25(0.15) 
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15.14(-0.76) 


15.78(-0.12) 


16.12(0.22) 






io«Cd 


17.2±0.3 


16.02(-1.18) 


16.63(-0.57) 
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"°Cd 


16.9±0.3 


15.74(-1.16) 


16.37(-0.53) 


16.74(-0.16) 


Es = 




"^Cd 


16.8±0.2 


15.62(-1.18) 


16.24(-0.56) 


16.61(-0.19) 






"^Cd 


16.7±0.4 


15.49(-1.21) 


16.11(-0.59) 


16.49(-0.21) 






"«Cd 


16.6±0.3 


15.35(-1.25) 


15.97(-0.63) 


16.37(-0.23) 



value of the nuclear matter incompressibility associated 
with this parameter set, namely Kqo = 201 MeV. Also 
in this case the results depend on the choice of the pair- 
ing force: when volume and mixed pairing interactions 
are adopted, the results tend to worsen. 

Whereas we can clearly confirm from the present re- 
sults that the value of the nuclear matter incompressibil- 
ity does play a key role in dictating the location of the 
ISGMR centroid energy, it is also true that the pairing 
interaction lowers the energy of the ISGMR to some ex- 
tent, typically few hundreds keV. This qualitative con- 
clusion is the same that was first found in Ref. Ill |. 



Thus, the pairing interaction cannot be neglected if one 
aims to reproduce not only the ISGMR centroid energies 
in Cd isotopes, but also, more generally, in other open- 



shell nuclei. It should also be noticed that the slope of 
the isotope dependence of the ISGMR centroid energy 
is rather well reproduced by all the three interactions, 
while the absolute values are much more sensitive to the 
choice of the Skyrme parameter set and of the pairing 
force. 

The various kinds of centroid energies are shown in 
Table II. The results that we report here, are obtained 
with the mixed pairing interaction for each Skyrme pa- 
rameter set. Again, the SkP interaction underestimates 
all the energies in Cd isotopes. With the mixed pair- 
ing force, SkM* can give an exc ellent des cription of the 
constrained energies Econ = a/ mi/m_i for all nuclei, 
while a slight underestimation of the centroid energies 
Ecen = 'Tii/mg is produced. Compared to the other two 
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FIG. 2: The calculated ISGMR strength distribution in los-necd are compared with the experimental data obtained at 
RCNP 7]. The SLy5 (solid line), SkM* (dashed line) and SkP (dotted line) forces are adopted in the calculations whose 
results are shown here, together with the mixed pairing interaction. 



interactions, in general, SLy5 can provide better results 
for the energies that are displayed in Table II. 

It is not completely clear, from the experimental point 
of view, whether the constrained, centroid or scaling 
energies are more suitable to be compared with the 
experimental data. However, from what we have just 
concluded, it can be stated that the reasonable values 
of nuclear incompressibility that can be extracted from 
the present Cd data are either the one of SkM* or SLy5, 
namely 217 or 230 MeV. This number is consistent with 



the one extracted from ^"^Pb, although on the lower 
side. 



B. Sn Isotopes 

Fig. 4 shows the calculated ISGMR strength distri- 
butions in the even-even ^^^~^^^Sn together with the 
experimental data taken from Refs. The results 

obtained with the Skyrme sets SLy5 (solid line), SkM* 
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FIG. 3; The calculated ISGMR centroid energies in the even-even los-uf'Qfj isotopes are compared with the experimental 
data obtained from Ref. The forces SLy5 (a), SkM* (b) and SkP (c) are adopted in the present calculations together 
with either the filling approximation, or the volume, or the surface, or the mixed pairing interactions, respectively. 



(dashed line) and SkP (dotted line), and with the mixed 
pairing interaction, are displayed. All three Skyrme in- 
teractions give a single peak around 16 MeV. The in- 
teraction dependence of the peak energy is qualitatively 
the same as in the case of Cd isotopes, i.e., the SkP re- 
sult is the lowest, the SkM* one lies in the middle and 
the SLy5 result is found at the highest energy - in agree- 
ment with the associated values of the incompressibility. 

The mass number dependence of the calculated IS- 
GMR centroid energies in the Sn isotopes is shown in 
Fig. 5, for the case of the SLy5, SkM* and SkP sets with 
the various pairing interactions. The results obtained 
without pairing, by using the filling approximation, are 
also shown in the same figure. It should be noticed that 
in some cases the experimental data are not consistent 
with each other. However, at the same time it is quite 
clear that the predicted centroid energies obtained using 
the SLy5 and SkP interactions are not in agreement with 
the experimental data: the SLy5 results, even when the 
pairing effects are taken into account, overestimate the 
experimental data while the SkP results underestimate 
them. In the case of the SkM* interaction, the predic- 
tions are much improved compared to the case of SLy5 
and SkP. The calculations with the volume and mixed 
pairing forces (and the filling approximation) give bet- 
ter predictions for ^^^'^^^Sn. The results provided by 



the volume pairing force in ii6>ii8,i20gj^ ^^^.g very close 
to the experimental data, while they slightly overesti- 
mate the experimental findings in ^^^'^^*Sn. 

There is some difference between the results presented 
here, and those already discussed in Ref. [ll[. Com- 
pared to Ref. [ll|, we have included here the con- 
tribution of the two-body spin-orbit interaction which 
provides an attractive effect and lowers the RPA and 
QRPA ISGMR energies. However, in order to be able 
to perform many systematic calculations, we have cho- 
sen in the current work the QRPA on top of HF-BCS 
instead of HFB. Pairing effects are slightly different in 
the two approaches, so the discussion of this Section on 
the results obtained with either volume, or surface, or 
mixed pairing forces, is not exactly the same as in Ref. 
[ll| . Despite these differences, we have to stress that 
the overall qualitative conclusion that the Sn data are 
rather consistent with the value of nuclear incompress- 
ibility associated with SkM* once the pairing is taken 
into account, is the same as in Ref. [TTj . 

This conclusion can be reinforced by the values re- 
ported in Table HI, that are the various calculated en- 
ergies obtained by means of the QRPA calculations per- 
formed with the mixed pairing interaction. The interac- 
tion SkP underestimates all these energies with respect 
to experiment, whereas SLy5 overestimates them. The 
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TABLE III: The calculated ISGMR constrained energies {Econ = -\/mi/m_i), centroid energies {Ecen = mi/mo), and 
scaling energies {Es — \J mzjmx) in even-even ^^^~^'^'*Sn isotopes are compared with the experimental data. The theoretical 
results are obtained in the interval between 10.5 and 20.5 MeV by using the SkP, Sk M*, and SLy5 parameter sets together 
with the mixed pairing interaction. The experimental data are taken from Refs. [sl. [gI. [soI. [sH . The values in parenthesis are 
the difference between the theoretical values and the experimental data. Units are MeV. 
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15.67(0.27) 


16.03(0.63) 
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15.74(-0.96) 


16.38(-0.32) 


16.77(0.07) 






"^Sn 


16.5±0.2 


15.62(-0.88) 


16.18(-0.32) 


16.66(0.16) 






"«Sn 


16.3±0.2 


15.53(-0.77) 


16.16(-0.14) 


16.54(0.24) 


Es = 




"«Sn 


16.3±0.1 


15.40(-0.9) 


16.05(-0.25) 


16.45(0.15) 






i2°Sn 


16.2±0.2 


15.29(-0.91) 


15.96(-0.24) 


16.35(0.15) 






i^^Sn 


15.9±0.2 


15.18(-0.71) 


15.86(-0.04) 


16.26(0.36) 






i^^Sn 


15.8±0.1 
14.96t°:l? 


15.07(-0.73) 


15.76(-0.04) 


16.18(0.38) 



interaction SkM* gives better results for both the con- 
strained energies, the centroid energies, and the scaling 
energies. 

Thus, the discrepancy between the values of the nu- 
clear incompressibility extracted either from Sn or Pb 
data still remains to some extent a puzzle; however, the 
pairing effects need to be taken into account and reduce 
this discrepancy to only « 5%. 



C. Pb Isotopes 

Measurements of the ISGMR strength distributions 
in the even-even nuclei 204,206p|^ have been done re- 
cently at RCNP, Osaka University [l^l- Several mea- 
surements of the ISGMR strength distribution in ^°*Pb 
had been performed already in the past, e.g., at KVI 
(Netherlands) and Texas A&M University (USA). One 
of the motivations to perform the new measurement on 
Pb isotopes at RCNP has been the study of a conjecture, 
that is, the possible appearance in the ISGMR energies 



of the so-called mutually-enhanced-magicity (MEM) ef- 
fect proposed by Lumiey and Zeldes in the context of 
the mass systematics [Sj, [11] . According to this con- 
jecture, the ISGMR energy in double closed shell nuclei 
might be higher if compared with the systematic values 
in its neighboring open shell isotopes [18|. Experimen- 
tally, the peak energies that have been obtained for the 
204-208pb isotopes [l^ are 13.98, 13.94, and 13.90 MeV, 
respectively. From these experimental results, we can 
not infer any kind of MEM effect in Pb isotopes. 

In Fig. 6 and Table IV we show our theoretical results 
for the ISGMR in Pb isotopes. The centroid energies 
calculated by using different Skyrme sets and pairing 
forces are displayed in Fig. 6. We can see that both 
SkM* and SkP underestimate the experimental data, al- 
though the results obtained with SkM* are better. The 
Skyrme force SLy5, having an associated value of in- 
compressibility Koo= 230 MeV, reproduce very well the 
centroid energies in Pb isotopes. The same conclusion 
can be drawn from the results shown in Table IV. It has 
been known for some time that this value of Koo can 
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TABLE IV: The calculated ISGMR constrained energies {Econ ~ \/mi/m-i), centroid energies (Ecen ~ mi/mo), and scaling 
energies {Es = -^ms/mi) in even-even 204,206,208p|^ compared with the experimental data. The theoretical results are 
obtained in the interval between 10 and 20 MeV by using the SkP, SkM*, and SLy5 parameter sets, together with the mixed 
pairing interaction. The experimental data are taken from Refs. [lol.l30t. The values in parenthesis are the difference between 
the theoretical values and the experimental data. Units are MeV. 









Exp. 


SKP 


SKM* 


SLy5 






204pi3 




12.84 


13.40 


13.80 




^mi/rri-i 


206 pj, 




12.79 


13.36 


13.76 






208 pj^ 




12.70 


13.29 


13.71 






204pij 


13.98 


12.88(-1.10) 


13.44(-0.54) 


13.84(-0.14) 


Ecen 


= mi/mo 


206 pj^ 


13.94 


12.83(-1.11) 


13.40(-0.54) 


13.81(-0.13) 






208 pb 


13.96±0.2 


12.74(-1.22) 


13.34(-0.62) 


13.77(-0.19) 










13.01 


13.56 


13.98 


Es = 




206 pj^ 




12.97 


13.54 


13.97 






208pt, 




12.88 


13.49 


13.93 



be extracted from the •^"^Pb data, if the density depen- 
dence of the force is the one that characterizes most of 
the recent Skyrme forces 3- Indeed, from the present 
results one can see that the effect of pairing in the Pb 
isotopes is rather small. Finally, there is no evidence for 
any MEM effect in our theoretical calculations which is 
also confirmed by the experimental results obtained at 
RCNP. 

IV. SUMMARY 

In summary, we have studied systematically the 
isoscalar giant monopole resonance (ISGMR) in Cd, 
Sn and Pb isotopes within the self-consistent Skyrme 
HF+BCS and quasi-particle random phase approxima- 
tion (QRPA). Three different Skyrme parameter sets are 
used in the present calculations, namely SLy5, SkM* 
and SkP: they are chosen since they are character- 
ized by different values of the nuclear incompressibility, 
Koo=230, 217 and 202 MeV, respectively. To study the 
role of the pairing correlations, we choose three types of 
pairing interactions, i.e., the so-called volume, surface, 
and mixed pairing forces. For the sake of comparison, 
we also produce RPA results (without pairing) within 
the filling approximation. 

The various kinds of centroid energies, and the de- 
tailed strength distributions of the ISGMR in Gd, Sn 
and Pb isotopes are compared with the available ex- 
perimental data. We have found that the pairing cor- 
relations always decrease the peak energies of the IS- 
GMR because of the attractive character of the particle- 
particle force in the 0+ channel. The typical size of this 
effect is several hundreds keV. 

From the present study, we find that the SkP in- 
teraction fails completely in reproducing the ISGMR 
strength distributions in all the studied isotopes due to 
its low value of nuclear matter incompressibility. The 
SLy5 parameter set (having an associated incompress- 
ibility of 230 MeV) gives a reasonable description of the 
ISGMR in Pb isotopes, whereas a better overall descrip- 
tion in the case of Gd and Sn isotopes is achieved by 
using the force SkM* (characterized by incompressibil- 



ity of 217 MeV). We have also found that the change of 
the ISGMR energies by pairing correlations in 204,206p|-, 
is quite small. The results for the Pb isotopes suggest 
that both theoretically and experimentally there is no 
evidence for the so-called MEM effect. 

Pairing helps in reducing the discrepancy between 
the values of the nuclear incompressibility extracted ei- 
ther from Pb or Sn data. We also found that Cd data 
do not introduce further problems as they seem to be 
rather consistent with Sn data. A small discrepancy of 
about 5% in Kqo between the conclusions drawn from 
Pb and Sn data remains and may deserve further in- 
vestigation. Since the size of this discrepancy depends 
on the pairing force and the model to treat its effects, 
one should better analyze whether we can constrain the 
attractive particle-particle matrix elements that appear 
in the QRPA calculations. Pairing forces are usually 
constrained only by means of ground-state properties, 
and at present it is unclear if volume or surface pairing 
interactions should be preferred. Along a different line, 
the present results may be interpreted by saying that 
the role of surface and surface-symmetry contributions 
to the nuclear incompressibility are still not precisely 
fixed. 



Acknowledgments 

We would like to thank M. Fujiwara, Y.-W. Lui and 
U. Garg for stimulating discussions and for providing 
us the experimental data prior to publication. This 
work is partially supported by the Japanese Ministry 
of Education, Gulture, Sports, Science and Technology 
by Grant-in-Aid for Scientific Research under the pro- 
gram number (G(2))20540277. L.G. acknowledges the 
support of the National Science Foundation of Ghina 
under Grant Nos. 10875150 and 11175216. The sup- 
port of the Italian Research Project "Many-body the- 
ory of nuclear systems and implications on the physics of 
neutron stars" (PRIN 2008) is also acknowledged. This 
research was supported in part by the Project of Knowl- 
edge Innovation Program (PKIP) of Ghinese Academy 



9 



2.5 



> 




E 

M— 

CO 

O 

CO 



0.0 



; (a) 


/% 


1 1 1 1 1 1 1 




•'W 


\ ■ Exp. 




* ■/ 

1 Jf 


•\ SLy5 : 


- 

1 


rl/ 


\& SkM* '■ 


t 

» 

t 

! 1. 


// 


'A SkP 


/x/i 
' n 




'•Ksl III 

^^dilD I i { 



12 16 20 24 28 
E(MeV) 



32 




16 20 24 
E(MeV) 



32 



2.5 



> 




1.5 



CO 



0.0 



: (c) 


1 1 1 1 1 1 1 

1 


1 . 1 


t\ 



8 12 16 20 24 28 
E(MeV) 



32 




16 20 24 
E(MeV) 



28 32 



2.5 

2.0 

$ 1.5 
E 

o 1.0 



> 




CO 



0.5 
0.0 



: (e) . 
• < 

1 


1 ' 1 ' 1 ' 1 ' 

'A 

• '\ 

• *\ 


^ ■'■ . 1 


• A 



> 





2.5 
2.0 
1.5 



8 12 16 20 24 28 
E(MeV) 



32 



"o 1.0 
CO 

0.5 
0.0 



; (f) 


■ 1 ■ 1 1 1 ■ 1 1 

' !\\ Sn 

•'4 

/ 'A 


T-' ' 

' . 1 

; // 

■ ■•'/ 

'5* 





8 12 16 20 24 
E(MeV) 



28 32 




8 12 16 20 24 28 32 
E(MeV) 



FIG. 4: The calculated ISGMR strength distributions in ^^^Sn are compared with the experimental data from Refs. 

F3 frll T^T CT ,.c: / — i;j i; — \ ci.ivi"* ^ i:„ 



,j — J- — ] 4-T~„ 1 i„ 



10 



16.5 



16.0 



>15.5 



15.0 
11 



1 — ' — I — ' — I — 
Sn isotopes 
SLy5 




■ Exp. RCNP • Exp. T&M 

-T— filling approx. —♦—volume pairing 
surface pairing — *•— mixed pairing 



112 114 116 118 120 122 124 
A 



16.5 



■^16.0 



>15.5 

15.0 
1 



Sn isotopea 
SkM* 




112 114 116 118 120 122 124 
A 



Sl6.0 




112 114 116 118 120 122 124 
A 



FIG. 5: The calculated ISGMR centroid energies in the even-even ^^''Sn isotopes are compared with the experimental 
data from Refs. The forces SLy5 (a), SkM* (b) and SkP (c) are adopted in the present calculations together with 

either the filling approximation, or the volume, or the surface, or the mixed pairing interactions, respectively. 



of Sciences, Grant No. KJCX2.YW.W10. 



[1] M.N. Harakeh and A.M. Van Der Woude, Giant Reso- 
nances: Fundamental High-Frequency Modes of Nuclear 
Excitations (Oxford University Press, Oxford, 2001). 

[2] J.R Blaizot, Phys. Rep. 64, 171 (1980). 

[3] S. Shlomo, V.M. Kolomietz, and G. Colo, Eur. Phys. J. 
A30, 23 (2006). 

[4] G. Colo, N. V. Giai, J. Meyer, K. Bennaceur, and P. 
Bonche, Phys. Rev. C70, 024307 (2004). 

[5] D. Vretenar, T. Niksic, and P. Ring, Phys. Rev. C68, 
024310 (2003). 

[6] J. Piekarewicz, Phys. Rev. C69, 041301 (R) (2004). 

[7] U. Garg, arXiv: 1101 .3125 

[8] T. Li et al, Phys. Rev. Lett. 99, 162503 (2007). 

[9] T. Li et al, Phys. Rev. C81, 034309 (2010). 
[10] M. Fujiwara, private communication. 
[11] J. Li, G. Colo, and J. Meng, Phys. Rev. C78, 064304 
(2008). 

[12] J. Piekarewicz, Phys. Rev. C76, 031301(R) (2007). 
[13] T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. C78, 
034318 (2008). 

[14] H. Sagawa, S. Yoshida, CM. Zeng, J.Z. Gu, and X.Z. 

Zhang, Phys. Rev. C76, 034327 (2007). 
[15] J. Piekarewicz and M. Centelles, Phys. Rev. C79, 



054311 (2009). 

[16] M. Centelles, S.K. Patra, X. Roca-Maza, B.K. Sharma, 

P.D. Stevenson, and X. Vinas, J. Phys. G: Nucl. Part. 

Phys. 37, 075107 (2010). 
[17] E. Khan, Phys. Rev. C80, 011307(R) (2009). 
[18] E. Khan, Phys. Rev. C80, 057302 (2009). 
[19] O. Bohigas, A.M. Lane, and J. Martorell, Phys. Rep. 

51, 267 (1979). 
[20] L. CapeUi, G. Colo, and J. Li, Phys. Rev. C79, 054329 

(2009). 

[21] V. Tselyaev, J. Speth, S. Krewald, E. Litvinova, S. 
Kamerdzhiev, N. Lyutorovich, A. Avdeenkov, and F. 
Griimmer, Phys. Rev. C79, 034309 (2009). 

[22] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Scha- 
eflfer, Nucl. Phys. A635, 231 (1998). 

[23] J. Bartel, P. Quentin, M. Brack, C. Guet, and H.-B. 
Hakansson, Nucl. Phys. A386, 79 (1982). 

[24] J. Dobaczewski, H. Flocard, and J. Treiner, Nucl. Phys. 
A422, 103 (1984). 

[25] D.J. Rowe, Nuclear collective motion (Methuen, Lon- 
don, 1970). 

[26] P. Ring and P. Schuck, The Nuclear Many Body Problem 
(Springer, Berlin, 1980). 



11 



14.0 



> 




.13.5 



13.0 



Pb isotopes SLy5 

■ Exp. RCNP • Exp. Lui. 

Exp. KVI —▼—filling approx. 

-♦— volume pairing surface pairing 
-*— mixed pairing 



204 



206 
A 



208 



14.0 



> 




E 



=13.5 



13.0 




14.0 



> 




'13.5 



13.0 ■ 




FIG. 6: The calculated ISGMR centroid energies in the even-even ~ Pb isotopes are compared with the experimental 
data obtained from Refs. [lol. [30l. [sH . The forces SLy5(a), SkM*(b) and SkP(c) are adopted in the present calculations 
together with either the filling approximation, or the volume, or the surface, or the mixed pairing interactions, respectively. 



[27] A.P. Severyukhin, Ch. Stoyanov, V.V. Voronov, and N. 

Van Giai, Phys. Rev. C66, 034304 (2002). 
[28] G. Colo, L. Cao, N. Van Giai, L. CapeUi, Comp. Phys. 

Comm. (submitted). 
[29] T. Sil, S. Shlomo, B.K. Agrawal, and P.-G. Reinhard, 

Phys. Rev. C73, 034316 (2006). 
[30] Y.-W. Lui, D.H. Youngblood, Y. Tokimoto, H.L. Clark, 

and B. John, Phys. Rev. C70, 014307 (2004). 



[31] D.H. Youngblood, Y.-W. Lui, H.L. Clark, B. John, 
Y. Tokimoto, and X. Chen, Phys. Rev. C69, 034315 
(2004). 

[32] D. Lunney, J.M. Pearson, and C. Thibauh, Rev. Mod. 

Phys. 75, 1021 (2003). 
[33] N. Zeldes, T.S. Dumitrescu, and H.S. Kohler, Nucl. 

Phys. A399, 11 (1983). 



